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Abstract. We give conditions under which the normalized marginal distri- 
bution of a semimartingale converges to a Gaussian limit law as time tends to 
zero. In particular, our result is applicable to solutions of stochastic differential 
equations with locally bounded and continuous coefficients. The limit theorems 
are subsequently extended to functional central limit theorems on the process 
level. We present two applications of the results in the field of mathematical 
04 ' finance: to the pricing of at-the-money digital options with short maturities 

and short time implied volatility skews. 
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1. Introduction 



Limit theorems for finite-dimensional stochastic processes as time goes to in- 
finity have been a classical object of study in probability theory and many results 
on the existence and uniqueness of invariant distributions, the convergence of the 
\ processes to the latter and the limiting behavior of the fluctuations around the 

limiting distributions have been obtained (see e.g. J2U], [22], [2E], [2H], [3D] and the 
references therein). More recently, small time asymptotics of finite-dimensional 
continuous time stochastic processes have attracted much attention. Apart from 
the theoretical interest, these have become of great importance in various ap- 
plied fields such as mathematical finance, where the increasingly high frequency 
of trades in financial markets requires pricing models behaving reasonably both 
on very short and on long time horizons. 

K 

In the works [2], [3], [I], [6], [13], [H], [21] and the references therein the 
authors study the behavior of the random variables E[/(A" 4o+ 5)|J r i ^] for small 
values of 5 > 0, where A" is a finite-dimensional (jump-) diffusion process, a Levy 
process or more generally a semimartingale, (J-" t x )i>o is the filtration it generates 
and the function / is taken from a space of suitable real-valued test functions. 
In [1], this program is carried out for general finite-dimensional semimartingales 
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and under appropriate continuity assumptions on the characteristics of X as well 
as smoothness assumptions on the function /, the almost sure limit 

(1.1) Mmr^El/fe)^] - f(X t0 )) 
is determined. 

Here, we are interested in small time Central Limit Theorems for finite di- 
mensional semimartingales; that is, instead of the almost sure limit (11. ip we are 
concerned with the limit 

(1.2) \im6-^{f(X 5 )-f(X )) 

in distribution. More precisely, we give sufficient conditions on the semimartin- 
gale X under which, for every suitable test function /, the limit (II. 2p exists and 
is given by a centered normal random variable (whose variance depends on the 
particular choice of the function /). The most closely related result in the litera- 
ture seems to be Theorem 2.5 of Doney and Mailer [11], which characterizes the 
Levy processes that satisfy a small time Central Limit Theorem. 

In addition to the just described Central Limit Theorems, we prove Functional 
Central Limit Theorems on the process level and give two applications of our 
results in the field of mathematical finance: to the pricing of digital options and 
the asymptotics of implied volatility skews. To outline the first of the two appli- 
cations, we recall that the price of a digital option with strike K and maturity 
5 on an underlying security with price process X in the presence of a constant 
interest rate r > is given by the formula 

(1.3) E[e~ rS 1 {Xs >k}} = e- r5 F(X s > K). 

In the limit 5 \ 0, that is for short maturities, this price tends to if K > Xq 
(out-of-the-money options) and to 1 if K < Xq (in-the-money options) as soon 
as X has right-continuous sample paths. The evaluation of the limit in the case 
K = Xq (at-the-money options) is however much trickier and, in general, the limit 
can take all values in the interval [0, 1] as we show below. However, if a Central 
Limit Theorem of the type described above holds for the semimartingale X, then 
the limit must be given by |. Moreover, in a special case we can bound the price 
in ( 11. 31) for any fixed value of 6 > from above and below by completely explicit 
functions tending to | in the limit 6 \ 0. By a well known relation between 
digital prices and implied volatility skews, we deduce bounds on the latter in 
certain models with stochastic interest rates. 

For the sake of a cleaner exposition, we first give the assumptions on the semi- 
martingale X and state our main results in the case of continuous trajectories. 

Assumption 1. Let T > 0, x G R m . Let X = (X t \ . . . , X^)J mT] be an R m - 
valued continuous semimartingale with canonical decomposition (see e.g. page 337 
in |21] j X — Xq = M + A, where M is a continuous local martingale, and A has 
locally finite variation. Assume that 



CENTRAL LIMIT THEOREMS 3 

(1) X = x a.s.; 

(2) there exists an a.s. positive stopping time ta such that a.s. 

A\ = [ bids, te [Q,t a ], 3 G {l,...,m}, 
Jo 

for an adapted process b; 

(3) there exists a random variable Cb, such that \bl\ < Cy, < oo for a.e. 
t G [0,ta] a.s., j G {1, . . . ,m}/ 

(4) there exists an a.s. positive stopping time tm such that the covariation is 
a.s. 

nt m 

{M*,M k ) t = / yVVfds, te[0,r M ], j,ke{l,...,m}, 

J ° 1=1 

for a progressive process a; 

(5) there exists a deterministic constant C a < oo, such that \a{ \ < C a for 
a.e. t G [0, t m ] a.s., j, k G {1, . . . , m}; 

(6) as £ \ ; o t — > L a.s., where L is a deterministic m x m-matrix; 

With this notation the Central Limit Theorem and the Functional Central 
Limit Theorem for continuous semimartingales read as follows. 

Theorem 2 (Central Limit Theorem). Let X satisfy Assumption^ Then for 
every f : W 71 — > MJ 1 such that there exists an open neighborhood U of Xq with 
f EC 2 (U,R n ), we have 

±=(f(X t )-f(x ))AN f ast\0, 

where Nf is a normal random vector with mean and covariance matrix 

V=(Df)(x )L(Df(x )L) T . 
Here, (Df)(xo) stands for the Jacobian of f at xq. 

Theorem 3 (Functional Central Limit Theorem). Let X satisfy Assumption^ 
Then for every f : R m — > IR n such that there exists an open neighborhood U of Xq 
with f G C 2 (U, M 71 ), the processes 

y/,u ;= f f{X«)-fM \ uG(Qjl)j 
V \ u / te[o,T] 
converge in law to a Brownian motion with variance- covariance matrix 

V=(Df)(x )L(Df(x )L) T . 

as u \ 0. 

We remark at this point that Assumption [1] is satisfied for weak solutions of 
stochastic differential equations (SDEs) under minimal regularity assumptions 
on the coefficients. 
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Remark 4. Let X be a weak solution of the m- dimensional SDE 

d 

dX{ = b 3 (t ) X t )dt + ^a jk (t ) X t )&B k t1 t>0, je{0,...,m}, 

k=l 

Xq = Xq O.S., 

where B is a standard d- dimensional Brownian motion, xq £ b : [0,T] x 
R m — > W m is uniformly bounded in a neighborhood of (0,xq) and a : [0, T] x 
R m — > R mxd is continuous in (Q,xq). Then, X satisfies Assumption [I] and, 
hence, Theorems [2] and [3] apply. 

We also note that, if X satisfies Assumption [1] and the matrix L is non-singular, 
then the price of an at-the-money digital option in (11.31) (that is, when K = xq) 
converges to \ in the limit 8 \, 0. This result can be significantly sharpened, 
when X is given by a weak solution of an SDE of the following type. 

Theorem 5. Suppose that the process X solves the stochastic differential equation 

(1.4) dX t = b(t,-)dt + a(t)dB u 

Xq = Xq, 

where b : [0, oo) x Q — > M. m is a bounded predictable process, a : [0, oo) — > M. mxm 
is a locally square integrable function taking values in the set of invertible matrices 
such that the smallest eigenvalue of cr(-) T er(-) is uniformly bounded away from 
and B is a standard m- dimensional Brownian motion. Then, the bounds 

(1.5) e /l(t) < ¥{X] > Xq 1 ) < e Mt \ t > 
apply. Here, the functions fi, fi are given by 



fl(t) = ~ i + W n, T ' Mog2 + 





1 &lll,ooHog2 



2 



where \\o~ 1 &||2,oo = su Pt, w W 1 (t)b(t,u)\2. Moreover, in the limit t \ ; the 
functions e^ 1 , e^ 2 admit the series expansions 



;i.6) e /i(t) = ^-y^f ik-^ikoo t^+o( t ), 



(1.7) ^«) = I + ^lk-U.oo^ + OCt). 

The rest of the paper is structured as follows. In sectional we give the proofs of 
Theorems [21 [5] and [3] in this order. In addition, we provide examples of continuous 
semimartingales, for which the limit in (11. 3p is not | and, therefore, the Central 
Limit Theorem (Theorem [5]) cannot hold with a non-degenerate Gaussian law in 
the limit. In section [3J, we state and prove extensions of Theorems [21 [5] and [3] to 
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semimartingales with jumps. Finally, in section HJ we explain the consequences of 
these results for the prices of at-the-money digital options with short maturities 
and the small time asymptotics of implied volatility skews. 



2. Continuous Semimartingales 
We start with the proof of Theorem [2j 

Proof of Theorem |5J Let / be as in the statement of the theorem and let Nf be 
an jV(0, V) random vector on some probability space (Q, A, P). We need to show 

<f(X t ) -f(x W 



(2.1) 



limE 

t\o 



9 



To this end, we fix a function g e Cb(IR n ,R), choose an open ball B such that 
B C U, and define the hitting time r := Tg c . Then with 

(2.2) 



t := r A t a A r M , 



we have 



E 



in — ^ — ) 



< 



E 



-EMNf)] 
f(X t )-f(x )^ ^ f(X tAT )-f(x ) ^ 

®r[9(Nf)] 



Vi 

E K /(XtAT ^ /(Xo) ) 

Hence in order to show (12. ip . it is sufficent to prove that the two summands 
in the latter upper bound tend to zero as t \ 0. Since the event {r = 0} 
has probability zero, the first summand converges to zero by the Dominated 
Convergence Theorem. Moreover, the convergence of the second summand to 
zero will follow, if we can show 



(2.3) 



f(X tAT ) - /(s ) d AT 

7t /5 



t\0. 



In order to prove (12. 3p . we first note that Doob's Integral Representation The- 
orem (see e.g. Theorem 18.12 on page 358 of [24]) in combination with part (j4]) 
of Assumption [1] implies the existence of an m-dimensional Brownian motion B 
(possibly on an extension of the primary probability space) such that a.s. 



(2.4) 



M 3 tAT 



m rtAr 

£ / <HB k s , te[0,T), j'e{l,..,m}. 
k=i Jo 

By part (J2]) of Assumption [1] and (12.41) we therefore have a.s. 

ptAT m rtAr 

(2.5) X{ AT = x + V s ds + J2 < k dB l te[0,T],je{ 
Jo k=1 Jo 



1 . . . , m}. 
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In addition, we recall that, by the Cramer- Wold Theorem, ( 12. 3j) holds iff for 

every s = (si, . . . , s n ) T G R™ 



(2.6) 



E 



fj(XtAr) - fj(x Q ) d ^~ 



Vt 



as t \ 0. To show this, we fix s = (si, . . . , s„) T G R™. Applying the local Ito 
formula (see e.g. Corollary 17.19 on page 341 of [24] ) in combination with (12.51) . 

we have with *5/ t = 04 )i<j,k<m '■= °tol f° r a ^ J 6 {!)•••) n \- 

/t/sr m rtAr Q r 

{£.f j ){X.)ds+Y, -f-{X s )af&B k s , t>0, 

where 



( £ *) = 5Etf^) + E^( u )' u (E U, se[0,r}. 



2 ^ T s <9a; fc 9a;z 

K,Z = 1 



fc=l 



Thus, we have for t > 0: 



(2.7) 



E fj(XtAr) — fj(Xp) _ J_ V~> 



tAr 



(£./j)Wdi 



flAT 



V 1=1 fc 7=1 ' 



By parts ([3]) and ((SJ) of Assumption [T] and the choice of B there exists a random 
variable C < oo a.s. such that sup ugB < C a.s. for s G [0, r], j G 

{1, . . . , n}. Thus, we have 



(2- 



-= Sj / {C s fi){X s ) ds < c Vt s j ^ °. * \ 

/l 3=1 J ° 3 = 1 



J= i 



a.s. 



Before examining the second term on the right-hand side of (12 .7p we observe 
that, for every t G [0, T], the random vector 



^Em=i^(^o)^^7 



is j\f(0, V) distributed. In particular, the distribution of Nt is independent of 



t, and N t = Nf for every t > 0. With ^ := |^(£o)£zfc we have for all 
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h G C h (R, 



E 



1 n m n 

h {-T f E a * E ^0 ] - Ep [/» (x; 

v i=i fc=i j=i 



fe=i 

n m 



< 



E 



+ 



1 

'^jE^E^ 

vr i=l fc=l 

n m 



J=L fc=l 
n 



n m 



.Jfc 

Mr 



j = l fc=l 
n m 



E 



k=l 
n m 



Tt TIL TL TTL 

(4 E s > E ^0 - (4= E a i E ) i{^<o 



<2||/i|| 0O P(r<t) ^0 
as t \ 0. Therefore, the random variables 



v j=i fc,z=i v i=i M=i 

converge in distribution to X]j=i s j^f as ^ \ 0. Next, we show that the difference 
between (12. 9p and the second term on the right-hand side of (12. 7p converges to 
zero in L 2 . By the Cauchy-Schwarz inequality and Ito's isometry we have 



E 



3=1 k,l 



al k - ^(x )L lk ) dB k s 



dx 



< nm 



E s ?E E 

j=i k,i=i 



j=l k,l=l 



< nm' 



E s ?E E 



1 

7 

t Jo 

t At 



rtAr 



(X s )a l s k -^(x )L lk )dB 



dxi y s> s dxi 



rtAr 



^(X>f-g M L,j-d, 



max 



( dfj (xW k 9fj 



(x )L 



ik 



L t se[0,tAr}\dXl 
3=1 k,l=l L - 

which indeed converges to zero as t \ by the Dominated Convergence The- 
orem. The just established L 2 convergence implies convergence in distribution. 
Summarizing, we have in the limit t \ 0: 

■y n rtAr 

-7iE s i / ( C sfj)(X s )ds ^0 a.s., 
Vt ~ Jo 

j=l k,l=l v u 3=1 

so that by Slutsky's theorem (12. 6p readily follows. □ 
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Remark 6. Recall that, if a family of probability measures satisfies a large de- 
viations principle (LDP) with a rate function I, then the validity of a CLT is 
related to the second derivative of I (see section 1.4 in [10J for a discussion in 
the case of Cramer's theorem). We now (heuristically) outline this connection in 
a very simple instance of our setup. Suppose that X satisfies a one- dimensional 
SDE (with zero drift for simplicity) 

X t = x + [ a(X s )dB s , t>0, 
Jo 

where a is bounded, bounded away from zero and Lipschitz continuous. Then, 
due to the time- change formalism for one- dimensional diffusions (see e.g. Theo- 
rem 8.5.1 on page 148 in [32],), for each 5 > 0, we can view the random variable 
X$ as the value of the diffusion 



Xf ] =x + V6 [ a(X^)dW s 
Jo 



at time 1, where W is the appropriate standard Brownian motion. Now, using the 
remark following Theorem 5.6.7 in [9] on page 214, and the contraction principle 
(see Theorem 4-2.1 on page 126 of ^\), we conclude that the random variables 
X t satisfy an LDP as £ \ with rate function 

1 f 1 f( s ) 2 

I{x + e) = - inf / — - ds. 

f(0)=x , 
f(l)=x +e 

Next, let £ be an antiderivative of the function 1/a. By the assumptions on a, 
the function £(/(•)) belongs to -f^ 1 ([0, 1]) if and only if the function f belongs to 
if 1 ([0, 1]). Hence, the latter infimum can be rewritten as 

inf / v(s) 2 ds. 

veHH[0,l]): Jo 
v(0)=S{xo), 
v(l)=T,(x +s) 

Due to Jensen's inequality, the infimum is reached when v is the affine function 
connecting S(xo) and S(xo + e). Plugging it in, we end up with 

(2.10) J(x + 5) = -[S(x + e)-S(x )] 2 = -^ -^J . 

That is, for e > small and fixed, we have the asymptotics 

(2.11) P(X t > x + e) ~ exp(-/(x + e)/t), 

where ~ stands for exponential equivalence. Now, pretend that we can apply the 
LDP ( 12. lip with a time- dependent e defined by e = z\fi, where z > 0. Since 
I(x ) = I'(x ) = 0, we have 

I(x + zVi) = \l"(xz)z 2 t + o(t), 
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and so 

P (^>,)^(__£_ +0(1) ). 

The Gaussian limit law is thus correctly identified by this heuristic argument (the 
case z < is similar). 

If a semimartingale X satisfies Assumption [H and the limit law in Theorem [2] 
is non-degenerate, we clearly have 

(2.12) lim F(X t > so) = \- 

We now give some examples where the value of this limit is not 1/2. 

Example 7. Let us consider the squared Brownian motion B 2 in one dimen- 
sion (no confusion with our superindex convention should arise). Then clearly 
lim t \oF(B 2 > 0) = 1, which does not contradict Theorem^ Indeed, the mar- 
tingale part in the canonical decomposition of B 2 is B 2 — t = 2 f B s dB s , which 
leads to 

(2.13) (B 2 - t) = 4 f B 2 ds -> 0, t \ a.s. 

Jo 

Since all items of Assumption^ are satisfied, Theorem^ tells us that -^B 2 con- 
verges in distribution to a degenerate normal random variable. 

Example 8. Denoting by $ the standard normal cumulative distribution func- 
tion, we see that for any p G (0, 1) and a standard Brownian motion B, the 
continuous process X t = B t + $ _1 (p)v / i satisfies F(X t > 0) = p for all t > 0. 
(Although not related to the present topic, we recall that the process B t = Wt + y/i 
occurs in Example 3.4 of Delbaen and Schachermayer |8J. They show that, when 
used as the price process of a financial security, X t (and also exp(X t ) ) allows for 
immediate arbitrage; the arbitrage disappears if proportional transaction costs are 
introduced [T9| Example 4.1] J 

The following example shows that each probability p e [0, 1) can even be 
realized by a continuous martingale. (Note also that that the non- continuous 
martingale t — P t , where P t is a Poisson process with parameter 1, satisfies 
lim tV) P(t-P t >0) = l.) 

Example 9. Consider the squared Bessel process of dimension 5 > 0, that is, 
the strong solution of the SDE 

dR s t = 2\fp~ldB t + 5dt 

with initial value Rq = 0. Then, the process Rf — 5t, t > is a martingale. We 
claim that 



for all p e [0, 1), there is a 5 G [0, oo) such that lim ¥(R s t - 5t > 0) = p. 
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The scaling property of squared Bessel processes (see section 1 in chapter XI of 
[55] ) shows 

lim F(R s t -5t>0) = P(i^ > 5). 

We show now that when one varies 5 in [0,oo) ; the right-hand side achieves all 
values p G [0, 1). For 5 > 0, the random variable Rf has the gamma distribution 
with shape parameter 5/2 and scale parameter 2 (see Corollary 1.4 in section 1 
of chapter XI in [33],). In particular, it has mean 5 and variance 25. We claim 
that 

(2.14) lim F(R( > 5) = 0. 

Let e > be arbitrary. By Chebyshev's inequality, we have 



F(R{ >5 + e) < 2 —, 



and so lim^o F{R\ > 5 + e) = 0. Therefore, recalling that, for any fixed < 5 < 
2, the density function of R^ is strictly decreasing, we have the estimates 

lim ¥(Ri >5) = lim P(5 + e > R{ > 5) 

<e\im 2- s l 2 Y{5/2)- 1 5 5 ' 2 ~ 1 e- 5 ' 2 

5 5 ' 2 - 1 _ e 
~ £ 8\o T{l + 5/2)/{5/2) ~ 2 

Thus, taking the limit e \ ; we end up with (12.141) . For 5 — > 00, the random 
variables (R{ - 5)/{25) 1 / 2 converge in distribution to a standard normal random 
variable [23j p. 340]. This implies 

lim F{R{ >5) = -. 

6— >oo 2 

It now follows from the Intermediate Value Theorem that, for every p 6 [0, \), 
we can find a 5 > such that 

lim F(Rf - 5t > 0) = p. 

(Note that for 5 = 0, we have R s t = 0, and so lim^o W(Rf - 8t > 0) = .) 
Finally, by considering the martingales 5t — Rf, we see that all values p G [0, 1) 
can be achieved. 

We now take a look at higher order terms beyond the limit in (12.121) . If 
X t = B t + bt is a one-dimensional Brownian motion with drift b G R, we have 
F(X t > x ) = \ + 0{t 1 / 2 ). Theorem El which we prove now, shows that this 
estimate persists for a larger class of Ito processes. 

Proof of Theorem O Fix a t > and make a change of probability measure 
according to the Girsanov Theorem, with the corresponding density being given 
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by 

Z -l ;= e -N t -l(N) t 



dP 



exp ( - J* atf-tys, •) dB s - ± jf ^(3)^(3, .)\\ da) . 



Under Q, the process X solves the equation 

(2.15) dX s = a(s) dB® 

on [0,t] with initial condition Xo = xq and where B 1 ^ is a standard Brownian 
motion under Q. Thus, Q(X| > Xq 1 ) = \ for all s G (0,t]. Moreover, 

(2.16) F(X} > X 1 ) = E« [Z, l m >x ol} ] • 

To obtain upper and lower bounds on the latter expression, we fix numbers 
p, q > 1 such that p~ l + g" 1 = 1 and apply Holder's inequality to deduce 

Q(Xl>Xji)=&[l {xl>xh} Zl'>Zr 1/p \ 

<E^l {Xtl>X(U ^] 1/p E«[V /P ] 1/9 - 
Taking the p-th power and rearranging, we get 

q(xl > xiy^[z;^]- p,q < ^[\xl>x h} Z t ] 

<Q(xl>X l ) 1/q E Q [Z?] 1/p , 

where the last upper bound follows again by Holder's inequality. This can be 
simplified to 



or 

(2.i7) (^f^[z;^- i r /q < nxi > x 1 ) < (i) 1/ V[^r i ] l/p - 

To estimate the bounds further, we note that 

z^/p-i = e - Wp +i)N t -± Wl +ir (N ) t . exp ^ + iy-( N ) t 

^ -{q/p+l)N t ~\( q /p+lf(N) t mavn (\(g 

and 



(2.18) < e -to/p+iW-|to/p+D a W« . exp + i)^^- 1 ^) 



z p-i = gCp-i^t-ICp-i) 2 ^* . exp (i(p _ i) p (7v) t ) 
(2.19) < e (p-i)^-5(p-i) 2 W* . exp {\{p - l)pt Wa-'bWlJ. 
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(Recall that we write ||o" _1 5||2,oo for sup t x |<t _1 &(£, a?) (2-) The first factors in ( 12. 18ft 
resp. (12.191) are P-martingales, since Novikov's condition is satisfied by our as- 
sumptions on b and a. Therefore, inserting these estimates into (I2.17p . we obtain 



p>i 



< nxi > x x ) 

'l\l/9 



< inf (-) " exp (1(^-1)^11^11^). 

p +q =1 

p>l 

It is easy to see that the lower bound is maximized by 



p = 1 + 



2 log 2 

whereas the upper bound is minimized by 



V 



2 log 2 



which together give ( 11. 5ft . Finally, the expansions given in the statement of the 
theorem can be computed by Taylor expansions of the explicit functions in the 
lower and upper bounds. □ 

We conclude this section with the proof of Theorem EJ 

Proof of Theorem [3j Let B be a Brownian motion with variance-covariance 
matrix V and let (ui)i & n be a sequence with elements in (0, 1) such that u\ \ 
as I — > oo. It is sufficient to verify the convergence of the finite-dimensional 
distributions 

(2.20) (Yl' Ul ,...,Y£ Ul )A(B tl ,...,BtJ, h,...,t w e[0,T], weN, 
and the tightness condition 

(2.21) lim lmT pf sup - > e) = 0, e > 0. 

<5\0 Z^oo V | 5 _ t [<5 ' / 

Indeed, by Theorem 1.3.2 in [35], condition ( I2.2ip implies the tightness of the 
laws of Y f ,u ' , I G M. Moreover, the convergence ( 12.20p allows to the identify the 
limit points with the law of B. 

First, we focus on (I2.20p . Fix t 1 , . . . , t w 6 [0, T] for some w G N; then by the 
Cramer- Wold theorem it suffices to show 

w n w n 

d=l j=l d=l j=l 
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for all s G M u,xri as I — > oo. Let r be defined as in (12. 2p . Arguing as in the proof 
of Theorem HI we see that it is enough to show 

w n w n 

d 



EE^sMEE**^ 

d=l j=l d=l j=l 

as I — y oo. However, this can be proven analogously to ( 12.61) . 



To show f)2.2ip . note that we may work with the stopped processes Y t /£ l , I G N. 
Indeed, since r, as defined in (12.21) . is a.s. positive, we have 

Hm pf sup \Y/' Ul - Y t f ™ I > e) = 0, e > 0. 

The triangle inequality thus shows that f 12 . 2 lj) is implied by 

(2.22) Jim lim" f( sup \(Y/£ l ) j - I >e)=0, e > 0, j G {1, . . . ,n}. 

By we get 
(2.23) 

P sup \(Y/^y-(Y t ^y\>e) <P sup — / \(£ r f j )(X r )\dr>~) 

\|«-t[<<5 / \\s-t\<8 y/Ul Ju t (sAT) l J 



+ p(sup f; r (tAT) ^-(X r )af dB* > £ -y^) 

\\s-t\<5 kv=1 Ju l (sAr) OXv * J 



According to ( 12.81) . we have 

»«j(Mt) 



— ► 0, £ > 0. 



Pf sup 4= ^ Ml dr > |) < p(c\/^ > 

We now investigate the second term on the right-hand side of (12.231) . After fixing 
5, j and I, we define the process 

In addition, we introduce the processes 

G\ := F l5+t - Fts, t G h := [0, 8], i G {0, . . . , [T/5J - 1}, 
and for i = [T/5\, 

Gf /5J := F [T/s]s+t - F [T/S]S , t G I [m := [0, T - [T/5\]. 

These are continuous local martingales and, thus, each of them can be represented 
as a time changed Brownian motion (see e.g. Theorem 18.4 on page 352 of [2"3]): 
G\ = WJ Gi ,. Moreover, the quadratic variation of G l can be bounded according 
to 

{G i ) t < 1 Clu l 5, tel u te{0,...,[T/5\}, 
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where < 7 < 00 only depends on m and the Jacobian of / on the ball B (see 
the paragraph preceeding (12.21) for the definition of the latter). Now, consider the 
event {sup| t _ s | <(5 \F t — F s \ > pp}. Clearly, on this event there exist s 0) to £ [0, T] 

such that I so — to I < $ and |F to — F so \ > Without loss of generality we may 
assume that < s < 5 < t < 25 (the other cases can be dealt with in the same 
manner). Then, either \F$ — F so \ > or \F to — F$\ > ^p. In the first case 

we get 

(2.24) ^ < \F S - F so \ < \F S0 - F \ + \F S - F \ < 2 sup \F r - F \. 

4 re [0,5] 



In the second case we have 
(2.25) -V 



< \F to - F 5 \ < sup |F,5 +r - F s \. 

re[0,5] 



These considerations show that on the event {sup| 4 _ s | <(5 \F t — F s \ > pp} there ex- 
ists an index i E {0, . . . |T/^J} such, that sup tg/ \G\\ > pp. Putting everything 
together we obtain 



P( sup 

^|s-t|<<5 



Uj(*AT) Qf 



<p(sup|G*| > 
< VP(sup|Gj|> 



£\/ui 



dx 

for at least one % 



> 



i=0 

< I j + I ) P 



sup 

0<r<7C2 Ui 5 
£ 2 



128 7 C2<5 



5^0 



> 0. 



Here, the last estimate follows from Bernstein's inequality (see e.g. Exercise 3.16 
on page 153 of [33J). We have established (I2.22p and, thus, the proof is finished. 

□ 



3. Semimartingales with Jumps 

This section is devoted to the extensions of Theorems [2], [3] and |5] to semimartin- 
gales with jumps. We start by stating the assumptions on the semimartingale 
X, which will replace Assumption [T] when jumps are present. 

Assumption 10. For a T > 0, let X = (X/, . . . , X t m ) t T g[0 T] be an R m -valued 
semimartingale with decomposition X = X c + J, such that 

(1) X c is a continuous semimartingale satisfying Assumption^ 
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(2) the process J is given by 

Jt= [ [ ip{s,z) (U{ds,dz) - fi{ds,dz)) + [ [ (p(s,z)U(ds,dz), 

JO J B x JO Ju m \Bi 

where Bi denotes the unit ball in W 71 , II is a Poisson random measure 
on [0, T] x M. m with compensator /i; the M. m -valued processes ip, if are 
predictable with respect to the filtration generated by II and 

rT 



E 



\ip(s,z)\ (i(ds,dz) 

IJQ JB 1 



< oo: 



(3) There exists an a.s. positive stopping time Tj such that 

E |n-/i|([0,i Atj] x Bi) =o(t 1/2 ) as t\0. 

We can now formulate the analogue of Theorem [2] in the case of semimartingales 
with jumps. 

Theorem 11 (Central Limit Theorem with jumps). Let X satisfy Assumption 
[TUl Then for every f : W m — > W 1 such that there exists an open neighborhood U 
ofx with f G C 2 (U,R n ), we have 

±=(f(X t ) - f(x )) 1> N f ast\0, 

where Nf is a normal random vector with mean and covariance matrix 

V=(Df)(x )L(Df(x )L) T . 

Proof. Let r > be such that the closed ball B r (xo) with radius r around xq is 
contained in U. Further, we denote by B r / 2 (xo) the closed ball with radius r/2 
around x and define the hitting time r := , 2 ( x . ) c - Finally, we introduce the 
stopping time 

(3.1) t :=t At a At m Atj 

and notice that r a.s. positive. Therefore, by the same argument as in the proof 
of Theorem [21 it suffices to show 

-L(f(X tAT )-f(x ))AN f as t\0. 

By Ito's formula in the form of Proposition 8.19 in [7], we have for all j G 
{l,...,n} andt G [0,T]: 

(3.2) 

HAt m ptAr 

axtAr) - fi(x ) = / (c s mx s ) d S +j2 -f-iXsK AB k s 

J ° k,l=l J ° X{ 



(3.3) + 



(3.4) 



ptAr p 

/ / (fj(X s _ + z)) - fj{X s J)) (n(ds, dz) - fi(ds, dz)) 
Jo Jb! 

ptAr p 

+ / / {f j {x a _ + <p(s,z))-f j (x a j))n{& s ,te)- 

Jo JM, m \B 1 
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Arguing as in the proof of Theorem [2} we see that the vector of terms on the 
right-hand side of (13. 2p . rescaled by converges in distribution to Nf as t \ 0. 

Thus, the theorem will follow if we can show that the terms (13.31) and ( 13. 41) . 
rescaled by 4=, converge to zero in probability as t \ 0. 

The term (13. 3j) . rescaled by can be decomposed into a sum 7) 1 + of the 
following two terms: 



i rtAr 



ViJo 



^ rtAr 



ViJo 
Then: 



(fj(X s ^+ip(s,z)) - fj(X^)) l{|^( S)2 )|< r/ 2} (n(ds,dz) -/i(ds,d^)), 
(fj( x s- +^(s,z)) - fj(X s ^)) l{\^( s ,z)\>r/2} (n(ds,dz) -//(ds,dz)). 



Eflr/fl < 2|l/l ^|°) IU E [|n - „| ([0, t A rj] x Bj 

which converges to zero as t \ by part (3) of Assumption [10J Moreover, since 
J a.s. has only finitely many jumps of absolute size greater than r/2 on every 
finite time interval, T t 2 converges to a.s. as t \ 0. 

Lastly, the term ( I3.4p . rescaled by converges to zero a.s. as t \ 0, since J 
a.s. has only finitely many jumps of absolute size greater than 1 on every finite 
time interval. □ 

As in the case of continuous semimartingales, the Central Limit Theorem can 
be strengthened to a Functional Central Limit Theorem, which in the presence 
of jumps reads as follows. 

Theorem 12 (Functional Central Limit Theorem with jumps). Let X satisfy 
Assumption [TUJ Then for every f : M. m — > M n such that there exists an open 
neighborhood U of xo with f £ C 2 (U,M. n ), the processes 



/ te[o,T] 

converge in law to a Brownian motion with variance- covariance matrix given by 

V = (Df)(x )L(Df(x )L) T 

as u \ 0. 

Proof. For each / and u as in the statement of the theorem, we write Q^ u for the 
law of the process Y* ,u on D([0, T], M. n ), the space of right-continuous functions 
on [0, T] having left limits; moreover, we denote by Q[ ,u the law of the continuous 
part of Y f ' u on C([0, T], R n ). We claim first that the family (Q / '") u6 ( ,i) is tight 
on D([0,T],R n ) if and only if the family (Q f c ' u ) u e{o,i) is tight on C([0, T], M n ) and, 
moreover, that the limit points of the two families are the same. 
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To prove the claim, it suffices to show that for every e > and j G {1, . . . , n}: 

(3.5) F( sup \{J t f, y\ > e) ->• as u \ 0, 

H£[0,T] ' 

where J^ ,u denotes the jump part of F^' M . Indeed, if this is the case, then every 
converging subsequence of (Q^' u )ue(o,i) i n ^([0> corresponds to a converg- 

ing subsequence of (Ql' u )«e(o,i) i n C([0,T],IR n ) and the limits of the two subse- 
quences have to coincide. Now, since the stopping time defined in ( 13. ip is a.s. 
positive, (13. 5p is implied by 

(3.6) p( sup |(j/*") j | > £, T > uT) as u \ 0. 

Furthermore, by Ito's formula in the form of Proposition 8.19 in |7J, we have on 
the event {r > uT}: 

(3-7) 

(j/'T'=4= r ! (/ J (X S _+V( S ^))-/,(X S _)) (U(ds,dz)-^ds,dz)) 

V U J0 JB 1 

(3.8) +JL T [ (f^X,. + V (s,z)) - f,(X s _))U(ds,dz). 

V u JO JR m \Bi 

As in the proof of Theorem [TTJ we decompose the integral on the right-hand side 
of (13.71) according to whether \if)(s, z)\ < r/2, or \if)(s, z)\ > r/2, and call the two 
resulting processes (J^"' 1 )-? and (J^'"' 2 )- 7 , respectively. Since the process (Jf' u ' i y 
is obtained by integrating a predictable process with respect to a compensated 
Poisson random measure, it is a square-integrable martingale. Thus, by Doob's 
maximal inequality, we have 



(sup \(j{^ 1 y\>e/2,r>uT) 



He[o,t] 

2 ^ 
< — 

Ey/U 



(uT)A- 

(/j-(X 8 _ + - /j(A s _))l { |^( Si2) |< r /2}n(ds,dz) 



where we wrote II for H — fi. Moreover, the same argument as in the proof of 
Theorem [Til shows that the latter upper bound tends to zero as u \ (by virtue 
of part (3) of Assumption [TO]) . Finally, since a.s. the process J^ ,u has finitely 
many jumps of size greater than r/2 on every finite time interval, the random 
variables sup te r or i |(J t 2 ) 3 '| converge to zero a.s. as u \ 0. In addition, by the 
same reasoning, the supremum over t G [0, T] of (13.81) tends to zero a.s. as u \ 
as well. Putting everything together, we end up with (I3.6p . finishing the proof of 
the claim. 



Lastly, one can proceed as in the proof of Theorem [3]to first show the tightness 
of the family ((5c' u )ue(o,i) on C([0, T], M") and to subsequently identify each of its 
limit points with the law of a Brownian motion with variance-covariance matrix 
V. In view of the claim above, this finishes the proof. □ 
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We conclude this section by stating and proving the analogue of Theorem [5] in 
the presence of jumps. 

Theorem 13. Suppose that the process X solves the stochastic differential equa- 
tion 

(3.9) dX t = b(t,-)dt + a(t)dB t + [ 4>(t,y)n(dt,dy), 

Xq = Xq, 

where b : [0, oo) x Q — y M m is a bounded predictable process with respect to the 
filtration of the standard m- dimensional Brownian motion B, a : [0, oo) — > ]R mxm 
is a locally square integrable function taking values in the set of invertible matrices 
such that the smallest eigenvalue of cx(-) T er(-) is uniformly bounded away from 
and ip is a predictable process with respect to the filtration of the Poisson random 
measure II . 

Suppose further that II is symmetric with respect to y (so that, in particular, 
its compensator vanishes) and that ipi(t,y) = —ipi(t,y) for allt > and y G M m 
with probability 1. Then, the bounds 

(3.10) e /lW < F(Xl > Xq 1 ) < e Mt \ t > 
of Theorem^ apply with the same functions f\, f 2 as there. 



Proof. We start by fixing a t > and changing the underlying probability measure 
P to an equivalent probability measure Q according to 

(3.11) ^ = exp ( - £ a(s)-'b(s, •) dB s - I jf \a(s)- l b(s, -)|| ds) . 

Then, in view of the independence of the continuous and the jump parts of X un- 
der P and the Girsanov Theorem, the process X solves the stochastic differential 
equation 

(3.12) dX s = a(s)dB®+ I ^(s,y)U(ds,dy), s e [0,t] 



with a standard Brownian motion B^ under Q and initial condition X = xq. 
Moreover, the random variables 

Ul l) := ( ta(s)dBf} and C/ t (2) := f [ ^(a, y) U(ds, dy) 
v Jo 71 Jo Jr™ 
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are independent under Q and their distributions r/^ and r]^ are symmetric. 
Hence, 

Q(X 4 1 >X 1 )=Q(f/i 1) + f/f ) >0) 

Q(ui 1] > -c) r/ 2) (dc) + / Q(ui 1] > -c) r/ 2) (dc) 

J-oo 
oo rO 

Q(ui 1] > -c) r/ 2) (dc) + / 1 - Q(U^ 1] > c) ?/ 2) (dc) 

J-oo 
1 

~ 2" 

From now on, one can follow the lines of the proof of Theorem [5] to finish the 
proof. □ 

4. Digital options and the implied volatility slope 

Suppose that the one-dimensional, positive process S models the price of a 
financial asset, and that P is the pricing measure. The riskless rate is r > 0. The 
holder of a digital call option with maturity T and strike K receives the payoff 
1{s t >k} at maturity. Digital options are peculiar in that the owner receives the 
full payoff as soon as they are only slightly in the money, as opposed to call 
options, say, which kick in gradually. By the risk-neutral pricing formula, the 
value of the digital call at time zero is 

D(K,T) := e- rT E[l {ST>K }} = e- rT F(S T > K). 

There is a considerable literature on short-maturity approximations for option 
prices. For OTM (out-of-the-money; So < K) or ITM (in-the-money; Sq > K) 
digitals, the first order approximation is clear: As soon as the underlying S is 
a.s. right-continuous at t = 0, the Dominated Convergence Theorem yields 

nmD(jr,r) = / ifs ° <A " < 0TM) 

t^o v ; H tfS >K (ITM). 

Finer information on the OTM decay (which trivially also covers the ITM be- 
havior) comes from small time large deviations principles for the underlying. 
E.g., see Forde and Jacquier [15] for the case of the Heston model and references 
about other diffusion processes. Our CLT-type results are useful in the ATM 
case So = K. As an immediate consequence of our limit theorems, we enunciate: 

Theorem 14. If the process S satisfies the assumptions of Theorem\ll\ (in par- 
ticular, if it satisfies those of Theorem [2] or Remark [4]), and the limit law is 
non- degenerate, then the limiting price of an at-the-money digital call is 1/2: 

(4.1) VmD{S ,T) = ±. 

This (intuitive) result captures virtually all diffusion-based models that have 
been considered (Black-Scholes, CEV, Heston, Stein-Stein, etc.). Although it 
seems to be new in its generality, in particular for jump processes, some special 
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cases can be inferred from the literature (see the comment at the end of this 
section) . 

The jump processes used in financial modelling are often Levy processes. It is 
clear that a compensated compound Poisson process will yield an (unrealistic) 
ATM digital price limit of either zero or one (see the remark before Example EJ). 
As for the infinite activity case, limit laws are not the appropriate way to get 
a result like ( 14. lj) . Doney and Mailer [IT] have determined all Levy processes 
that admit a short-time CLT, with a criterion involving the tail of the Levy 
measure. While there do exist infinity activity Levy processes that satisfy a CLT 
[TTJ Remark 9], the Levy processes that have been considered in mathematical 
finance are typically not of this kind. For instance, it is easy to see from the 
characteristic function that the variance gamma process [27] does not admit any 
non-degenerate limit law for t — y 0, for any normalization. We will discuss these 
issues further in the more application-oriented companion paper 



Finally, we discuss the implied volatility skew. Suppose that the underlying S 
generates the call price surface C(K,T): 

C{K, T) = e~ rT E[(S T - K) + ], K > 0, T > 0. 

Then the implied volatility (see e.g. [25J) for strike K and maturity T is the 
volatility <Ti mp (K,T) that makes the Black-Scholes call price equal to C(K,T): 

C BS (K,a imp ,T) = C(K,T). 

The map K — y <Ji mp (K,T) is called the volatility smile for maturity T. It is also 
called the volatility skew, because it is often monotone instead of smile-shaped, 
but we will reserve the term skew for the derivative dRO~i mp (K,T). If C(K,T) is 
smooth in K, it equals (we omit arguments) 

o 9rCbs — 9rC 

OrO 'imp . 

Co-G-BS 

Under mild assumptions (e.g., if the law of St is absolutely continuous), we have 

(4.2) d K C = -e- rT F{S T > K) = -D(K,T), 

from which we deduce the (well-known) connection between the volatility skew 
and the price of a digital call (see e.g. [T7]): 

_ D(K,T) + d K C BS 

ORO-imp ~ T-^ • 

Co-G-BS 

Inserting the explicit Black-Scholes vega and digital price (see e.g. [31]), we obtain 

-D(K,T) + <5>(-a imp VT/2) 



9rO~ imp 



KVfn(a imp Vf/2) 



with $ and n denoting the standard normal cdf and density, respectively. For 
T — y 0, we have a imp \/T = o(l) under the following mild assumptions [3U 
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Proposition 4.1]: 

(4.3) (Sq - K) + < C(K, T) < S (no arbitrage bounds), 

(4.4) limC{K,T) = (S -K)+, 

T-»0 

(4.5) T i-> C(K, T) is non-decreasing. 
Therefore, 

(4-6) d K a imp r,^^--D(K,T)-^^^ + 0((a imp VT) 5 )), T -+ 0. 



We see that the small time behavior of the skew is related to that of the digital 
price. At the money, the latter will typically tend to 1/2 (see Theorem [T4"|) . and 
so higher order estimates are needed to get the first order asymptotics of the 
ATM skew dKO'i rap \K=s - To this end, we apply our Theorem [5j and compare our 
findings with the standard model free slope bounds [16j page 36] 



(4.7) - -^(1 - m))e~ rT+d " /2 < ^ < £^Hd 2 )e~^, 



where 



\og(S Q /K) + (r + la? mp )T 



0"imp\M 



d 2 = d 1 - a imp VT. 

Such bounds can give guidance on model choice; recall that the market slope 
seems to grow like T -1 / 2 for short maturities [I]. Note that the following result 
accomodates stochastic interest rates, and recall that we assume in this section 
that the dimension is m — 1 . Under stochastic interest rates, the digital call 
price is 

(4.8) D(K,T) = E[e-£r(s)^ 1{ST>K} y 

To calculate the implied volatility, a deterministic rate r has to be chosen (e.g., by 
e~ rT = E[exp(— J Q T r(s) ds)]). This choice is irrelevant for Theorem ITS] though. 

Theorem 15. Assume that the price process satisfies the SDE 

dS t /S t = r(t) dt + a(t) dB t 

with the stochastic short rate process [r(t)) t >o, o,nd that the log-price X = log S, 
whose drift is b(t) = r(t) — |<r 2 (t), satisfies the assumptions of Theorem [5j 
Assume further that d K C(K,T) = -D(K,T) holds (cf. (IP]) ) and that fOl) 
and (14. 5p are satisfied. Then we have the ATM slope bounds 



d K a^ p \ K=So > -^L -CVf - a -^J- + 0(T) + 0((a imp Vfr 
Kyi \ Ay An 



d K a imv \ K=So < ^0= (cVf - + 0(T) + 0((a imp Vff) j , 
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where 

Proof. According to (I4.8p . the ATM digital price equals 

D(S ,T)=E[e-^ T ^l {XT>xo} ]. 
The discount factor is 1 + O(T), so we can apply Theorem [5] to conclude 

~ - CVf + O(T) < D(S , T )<\ + CVT + 0(T). 

Now, the result follows from ( I4.6p . Note that cr imp \/T = o(l) by [MJ Proposi- 
tion 4.1], since we assume f !4.4j) and (14. 5p . and (14. 3[) is satisfied in our setup. □ 

The bounds in Theorem [15] are asymptotically stronger than the general esti- 
mate (14. 7p . which is of order 0(T~ 1//2 ), since ai mp y/T = o(l). If the Berestycki- 
Busca-Florent formula [5 J holds, then implied volatility tends to a constant. 
Therefore, our bounds are considerably stronger than ( 14. 7ft in this case, namely of 
order 0(1). The models covered by Theorem [151 thus do not match the empirical 
slope behavior T _1//2 , similarly to stochastic volatility models [26], whose slope 
also behaves like 0(1). 

To conclude our discussion of ATM digitals and the implied volatility skew, 
note that, for some diffusion processes, the result in Theorem [TH is implicitly in 
the literature. To wit, by (14. 6p . a non-exploding ATM slope requires a limit price 
of 1/2 of the digital. See Durrleman [12] page 59] for a general expression for the 
implied volatility slope that shows that it does not explode, e.g., in the Heston 
model. 

References 

[1] E. Alos, J. A. Leon, M. Pontier, and J. Vives, A Hull and White formula for 
a general stochastic volatility jump- diffusion model with applications to the study of the 
short-time behavior of the implied volatility, J. Appl. Math. Stoch. Anal., (2008), pp. Art. 
ID 359142, 17. 

[2] E. Alos, J. A. Leon, AND J. Vives, On the short-time behavior of the implied volatility 
for jump-diffusion models with stochastic volatility, Finance Stoch., 11 (2007), pp. 571-589. 

[3] E. Benhamou, E. Gobet, and M. Miri, Smart expansion and fast calibration for jump 
diffusions, Finance Stoch., 13 (2009), pp. 563-589. 

[4] A. Bentata and R. Cont, Short-time asymptotics for marginal distributions of semi- 
martingales, arXiv:1202.1302, (2012). 

[5] H. Berestycki, J. Busca, and I. Florent, Asymptotics and calibration of local volatil- 
ity models, Quant. Finance, 2 (2002), pp. 61-69. Special issue on volatility modelling. 

[6] H. Berestycki, J. Busca, and I. Florent, Computing the implied volatility in sto- 
chastic volatility models, Comm. Pure Appl. Math., 57 (2004), pp. 1352-1373. 

[7] R. Cont and P. Tankov, Financial modelling with jump processes, Chapman & 
Hall/CRC Financial Mathematics Series, Chapman & Hall/CRC, Boca Raton, FL, 2004. 

[8] F. Delbaen and W. Schachermayer, The existence of absolutely continuous local 
martingale measures, Ann. Appl. Probab., 5 (1995), pp. 926-945. 



CENTRAL LIMIT THEOREMS 



23 



[9] A. Dembo and O. Zeitouni, Large deviations techniques and applications, vol. 38 of 
Stochastic Modelling and Applied Probability, Springer- Verlag, Berlin, 2010. Corrected 
reprint of the second (1998) edition. 

[10] F. den Hollander, Large deviations, vol. 14 of Fields Institute Monographs, American 
Mathematical Society, Providence, RI, 2000. 

[11] R. A. Doney and R. A. Maller, Stability and attraction to normality for Levy processes 
at zero and at infinity, J. Theoret. Probab., 15 (2002), pp. 751-792. 

[12] V. Durrleman, From implied to spot volatilities, PhD thesis, Princeton University, 2004. 

[13] J. Feng, M. Forde, and J. -P. Fouque, Short-maturity asymptotics for a fast mean- 
reverting Heston stochastic volatility model, SIAM J. Financial Math., 1 (2010), pp. 126- 
141. 

[14] J. E. Figueroa-Lopez and C. Houdre, Small-time expansions for the transition dis- 
tributions of Levy processes, Stochastic Process. Appl., 119 (2009), pp. 3862-3889. 

[15] M. Forde and A. Jacquier, Small-time asymptotics for implied volatility under the 
Heston model, Int. J. Theor. Appl. Finance, 12 (2009), pp. 861-876. 

[16] J. -P. Fouque, G. Papanicolaou, and K. R. Sircar, Derivatives in financial markets 
with stochastic volatility, Cambridge University Press, Cambridge, 2000. 

[17] J. Gatheral, The Volatility Surface, A Practitioner's Guide, Wiley, 2006. 

[18] S. Gerhold and M. Shkolnikov, Short maturity implied volatility skew and digital 
option prices. Working paper, 2012. 

[19] P. Guasoni, No arbitrage under transaction costs, with fractional Brownian motion and 
beyond, Math. Finance, 16 (2006), pp. 569-582. 

[20] R. Z. HASMINSKH, Stochastic stability of differential equations, vol. 7 of Monographs 
and Textbooks on Mechanics of Solids and Fluids: Mechanics and Analysis, Sijthoff & 
Noordhoff, Alphen aan den Rijn, 1980. Translated from the Russian by D. Louvish. 

[21] J. Jacod, Asymptotic properties of power variations of Levy processes, ESAIM Probab. 
Stat., 11 (2007), pp. 173-196. 

[22] J. Jacod and A. N. Shiryaev, Limit theorems for stochastic processes, vol. 288 of 
Grundlehren der Mathematischen Wissenschaftcn [Fundamental Principles of Mathemat- 
ical Sciences], Springer- Verlag, Berlin, second ed., 2003. 

[23] N. L. Johnson, S. Kotz, and N. Balakrishnan, Continuous univariate distributions. 
Vol. 1, Wiley Series in Probability and Mathematical Statistics: Applied Probability and 
Statistics, John Wiley & Sons Inc., New York, second ed., 1994. A Wiley-Interscience 
Publication. 

[24] O. Kallenberg, Foundations of modern probability, Probability and its Applications 

(New York), Springer- Verlag, New York, second ed., 2002. 
[25] R. W. Lee, Implied volatility: statics, dynamics, and probabilistic interpretation, in Re- 
cent advances in applied probability, Springer, New York, 2005, pp. 241-268. 
[26] A. L. Lewis, Option valuation under stochastic volatility, Finance Press, Newport Beach, 

CA, 2000. With Mathematica code. 
[27] D. Madan, P. Carr, and E. Chang, The variance gamma process and option pricing, 

European Finance Review, 2 (1998), pp. 79-105. 
[28] S. P. Meyn and R. L. Tweedie, Stability of Markovian processes. I. Criteria for 

discrete-time chains, Adv. in Appl. Probab., 24 (1992), pp. 542-574. 
[29] , Stability of Markovian processes. II. Continuous-time processes and sampled chains, 

Adv. in Appl. Probab., 25 (1993), pp. 487-517. 
[30] , Stability of Markovian processes. III. Foster-Lyapunov criteria for continuous-time 

processes, Adv. in Appl. Probab., 25 (1993), pp. 518-548. 
[31] M. MUSIELA AND M. RUTKOWSKI, Martingale methods in financial modelling, vol. 36 of 

Stochastic Modelling and Applied Probability, Springer- Verlag, Berlin, second ed., 2005. 
[32] B. Oksendal, Stochastic differential equations: An introduction with applications, 

Springer- Verlag, Berlin, sixth ed., 2010. 



24 



S. GERHOLD, M. KLEINERT, P. PORKERT, AND M. SHKOLNIKOV 



[33] D. Revuz and M. Yor, Continuous martingales and Brownian motion, vol. 293 of 
Grundlehren der Mathematischen Wissenschaftcn [Fundamental Principles of Mathemat- 
ical Sciences], Springer- Verlag, Berlin, third ed., 1999. 

[34] M. Roper and M. Rutkowski, On the relationship between the call price surface and 
the implied volatility surface close to expiry, Int. J. Theor. Appl. Finance, 12 (2009), 
pp. 427-441. 

[35] D. W. Stroock and S. R. S. Varadhan, Multidimensional diffusion processes, Classics 
in Mathematics, Springer- Verlag, Berlin, 2006. Reprint of the 1997 edition. 

Vienna University of Technology, Wiedner Hauptstr. 8/105-1, A-1040 Austria 
E-mail address: sgerhold@f am.tuwien.ac.at 

Vienna University of Technology and arithmetica versicherungs- und finanz- 
mathematische beratungs-gmbh 

E-mail address: kleinert.max@gmail.com 

Christian Doppler Laboratory for Portfolio Risk Management, Vienna Uni- 
versity of Technology, Wiedner Hauptstr. 8/105-1, A-1040 Austria 
E-mail address: piet .porkert@f am. tuwien. ac . at 

Department of Statistics, University of California, Berkeley, CA 94720-3860 
E-mail address: mshkolni@gmail.com 



